Generalization of k ■ p theory for periodic perturbations 
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We extend standard k • p theory to take into account periodic perturbations which are rapidly 
oscillating with a wavelength of a few lattice constants. Our general formalism allows us to explicitly 
consider the Bragg reflections due to the perturbation-induced periodicity. As an example we 
calculate the effective masses in the lowest two conduction bands of spontaneously ordered GaInP2 
as a function of the degree of ordering. Comparison of our results for the lowest conduction band 
to available experimental data and to first principle calculations shows good agreement. 
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I. INTRODUCTION 



For many years k • p theorjfli has been very successful 
in describing a wide variety of crystal band structures, an 
important merit being that it allows one to derive simple, 
analytical formulas which capture the essential physics. 
In the presence of perturbing potentials it becomes 
the envelope function approximation (EFA),cl which has 
been successfully applied to such different problems as 
impuritieso and semiconductor heterostructures.u How- 
ever, the EFA requires the perturbation potential Vi(r) 
to be slowly varying on the length scale of the lattice 
constant, i.e., the nonzero Fourier components of the po- 
tential Vi (r) must be restricted to wave vectors which are 
small compared_to the dimensions of the Brillouin zone. 
Some authoraj'lD reported problems with the EFA for 
systems where this requirement is not fulfilled, notably 
in artificial and natural short-period superlattices such 
as spontaneously ordered GaInP2. These superlattices 
can be viewed as systems with a periodic perturbation, 
where the smallest Fourier components of the perturba- 
tion are comparable to the dimensions of the Brillouin 
zone of the unperturbed problem. These Fourier compo- 
nents result in interactions between states with very dif- 
ferent wave vectors, which are crucial for the properties 
of these short-period superlattices. However, the inter- 
actions canriot be accurately described within standard 
k ■ p theory.B In this paper we present a general method 
which allows us to treat perturbing potentials which are 
rapidly oscillating but commensurate to the periodicity 
of the potential of the unperturbed problem. To illus- 
trate the method, we apply it to the natural superlattice 
of spontaneously ordered GaInP2 . 

The Ga2;Ini_3;P alloy for x k, 0.51 can be lattice 
matched grown on a GaAs(OOI) substrate. Under proper 
growth conditions, long-range order of the CuPt type 
is observed.Hi3 This type of ordering is characterized by 
(III) layers alternatingly rich in Ga or In. In the ordered 
material the symmetry of the lattice is reduced from Td 
to Cat,, and the Brillouin zone becomes smaller than 
the zinc-blende Brillouin zone, which leads to a back- 
folding of states. The change in the crystal potential 



induced through the ordering is a short-period poten- 
tial that mixes electronic states of an "averaged" zinc- 
blende structure. In particular, the interactions between 
r and L states lead to energy shifts of the band-edge 
states in the ordered alloy, which cause band-gap reduc- 
tion and valence-band spliltiag. These effects hav 
investigated theoreticallyllj'y and experimentally. 
In addition to the changes in the energies ,|-the effective 
masses are also altered. Raikh and TsiperliJ calculated 
the conduction-band effective mass of ordered GaInP2 
using a two-band model, which accounts only for the mix- 
ing of conduction-band P and L states. They found that 
the effective mass parallel to the ordering direction m|| 
and the effective mass perpendicular to the ordering di- 
rection increase with increasing ordering, and that 
m|| is larger than m^. This model, however, does not 
take into account the change of the interaction between 
conduction- and valence-band states due to the band- 
gap reduction and valence-band splitting. These changes 
were investigated by Zhang and MascarenhasEj with an 
eight-band k • p model, which included zinc-blende P 
states from both the conduction and valence band. They 
find my and mj^ to decrease with increasing ordering, 
with TO II being larger than to^. A third investigation was 
done by Franceschetti, Wei, and ZungerJJ who performed 
first-principle calculations using the local-density approx- 
imation. They find TO|| to increase, whereas toj_ de- 
creases with increasing ordering. This result agrees qual- 
itatively with the only measurment that investigated t 
anisotropy of the effective masses in ordered GaInP2. 
In Ref. 1^ the conclusion was drawn that the conduction- 
band effective masses in ordered GaInP2 depend on a 
"delicate balance" of P-L mixing and increased interac- 
tion between conduction and valence band. However, the 
P-L mixing and the increase of the interaction between 
conduction and valence bands have a common source in 
the ordering-induced P~L interactions. The conduction 
band masses in partially ordered GaInP2 present an ex- 
cellent test for our theory, which should be able to de- 
scribe that "delicate balance." 

This paper is organized as follows. In Sec. || we de- 
rive a general scheme applicable to periodic perturbations 
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within k • p theory. In Sec. [11 we use this theory to de- 



rive a model for the conduction band of spontaneously 
ordered GaInP2. The results are discussed in Sec. ^M- A 
summary and a short outlook are presented in SecT|V[. 



II. GENERAL THEORY 

We consider a system with the one-particle Hamilto- 
nian 



(1) 



where p is the momentum operator, m the free-electron 
mass, and Vb(r) the periodic potential of the crystal. In 
the case of the GaInP2 alloy, Vb(r) is an averaged po- 
tential of the disordered material. The eigenfunctions 
of the Hamiltonian Hq are Bloch functions ipnk{^) — 
e''' ''u„k(r) with eigenvalues £„(k) 



(2) 



Here k is a wave vector in the first Brillouin zone (BZq), 
which corresponds to the periodicity of Vb(r). 
In conventional k ■ p theory, basis functions 



_ jk r iko-r 



Wnko(r), 



(3) 



are introduced, where k belongs to BZq. As shown by 
Luttinger and Kohn,El the functions x^k provide a com- 
plete and orthonormal basis set for any ko being an el- 
ement of BZq. Thus any eigenfunction of Hq can be 
expanded in terms of 

^n"ir)=J2 I d'k'c^„,A^')xX(.^), (4) 

BZo 

which yields the well-known k • p equation for expansion 
coefficients c„//„(k) and energy £„//(ko -I- k) 

J2 i?^°,(k)c„„„,(k) = e„.(ko + k)c„.„(k), (5) 

with 



HZ'i^) = { £n(ko) + d„„, + -k . p::^;,. (6) 



The momentum matrix element p„° / is defined as 



no 



d3re~^''°'-<koPe"'°"Wn'ko- 



(7) 



Here the integration extends over the unit cell with vol- 
ume rio- Equation (^) is diagonal with respect to k due 
to the periodicity of Vb(r), i.e., k is a good quantum 
number and the eigcnstates 0„ can be written as a su- 
perposition of basis functions x^w for different bands n 
but the same wave vector k. Because of the complete- 
ness of the basis functions x!^k('") arbitrary potential 
Vi (r) can be taken into account in Eq. (||) , even if Vi does 
not have the periodicity of Vq. However, this results in 
eigenstates 0„(r) which are superpositions of functions 
Xnk f*-"^' different k. In practical calculations one is usu- 
ally restricted to wave vectors k -I- ko close to ko, since 
the k ■ p term is treated as a perturbation. Therefore, 
this method will fail if the additional potential Vi(r) is 
not smooth. 

Here we will generalize the above approach to perturb- 
ing potentials Vi (r) which are rapidly oscillating but pe- 
riodic and commensurate to the periodicity of Vq. In 
the case a£ GaInP2, Vi(r) corresponds to the ordering 
potential.tS 

The total potential Vq -f Vi is characterized by a larger 
unit cell, and hence a smaller Brillouin zone (BZi). Due 
to commensurability, the reciprocal-lattice vectors asso- 
ciated with the original potential Vq can be expressed 
as an integer linear combination of the reciprocal-lattice 
vectors associated with the perturbing potential. In par- 
ticular we have several wave vectors in the larger Bril- 
louin zone BZo which become reciprocal-lattice vectors of 
the perturbed problem, i.e., these wave vectors become 
equivalent to the T point. This set of wave vectors will 
be called {/C}. For ordered GaInP2 this consists of the 
F and L point of BZq. 

The main difference from standard k • p theory is that 
eigenfunctions of -ffo belonging to wave vectors of the set 
{K} are used to form basis functions of the form 



X:fK(r) = e'^->„^:(r) 



iK-i 



(8) 



Any function having the periodicity of the perturbed sys- 
tem can be expanded in terms of the functions ^n7c(r) as 
the set {/C} is folded onto the F point of BZi. Therefore, 
when K is taken from BZi the functions X5w(r) form a 
complete and orthonormal basis (cf. also Refs. y|and [L?! ). 
The eigenfunctions oiH = Hq -\- Vi are expanded in terms 

<i>„„(r)=^ J d'K'cSn'i^')Xn'K'(.r). (9) 

This yields the Schrodinger equation for the expansion 
coefficients C^,„(K), 



— K P'^ 

m 



(10) 



where 



2 



^ ~lh~ J "«A:Pe Un'K (11) 

and 

On the left-hand side of Eq. ( [l0|) t he part diagonal in /C is identical to H^^, defined in Eq. (||), because we have 
-f*nn' ~ Pm'- [Note that in Eq. (O) the larger normalizing volume fli is compensated for by the larger integration 
volume.] The only effect of the perturbing potential Vi are the coupling matrix elements , and Eq. (|l^) can be 

written as 

^ff^„,(K)C^,„,(K) + J2 V,trC^:^,{K) = i?„„(K)C^,„(K). (13) 

n' n'JC' 



Equation (|lC|) represents our generalization of the stan- 
dard k • p equation Like^Eq. (|^), this new equation 
depends only explicitly on KJl3 

As noted above, the set {K} is folded onto the F point 
of BZq. Therefore, an expansion in terms of the basis 
functions (^) is appropriate for the description of states 
near the center of the Brillouin zone of the perturbed 
system. If one is interested in states near Ko ^ 0, one 
has to replace Khy K + Kg in Eq. (||) . 

It is straightforwar|d,to include spin-orbit coupling and 
strain- induced effectalj in Eq. (p^). Likewise, we can eas- 
ily add a slowly varying perturbation, like the potential 
of an impurity, as ansatz (w does not require periodicity. 



III. APPLICATION TO PARTIALLY ORDERED 
GalnPa 

In this section we apply the above theory to the 
conduction-band effective masses at the F point of par- 
tially ordered GaInP2. The unperturbed potential Vq 
corresponds to the disordered material. The perturba- 
tion Vi represents the ordering potential which is defined 
in Ref. as the difference between the potentials of the 
ordered and disordered material. In principle there are 
four equivalent variants of CuPt ordering for GaInP2, 
corresponding to the four (111) directions. Due to sub- 
strate effects onlij two of them are observed in exper- 
iments, however.H Since we consider the bulk system, 
the domain structure is not relevant to our calculations. 
Hence we choose the [111] direction to be parallel to the 
ordering direction. 

Due to the periodicity of the ordering potential Vi, 
zone center states of ordered GaInP2 are derived from 
zinc-blende F- and L-point states of the disordered sys- 
tem. Thus we choose {F, L} as the wave-vector set {K}. 
We restrict ourselves to a seven-band model, containing 
the zinc-blende Fic, Fsv, Lie, aad Lav states (nomen- 
clature according to Koster et allES). Spin-orbit interac- 
tion is neglected, as it has only a minor influence on the 
conduction-band effective masses. Due to time-reversal 
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symmetry wave functions from F and L points can be 
chosen to be real. With this phase convention all mo- 
mentum and potential matrix elements can be defined as 
real quantities. 

A. k ■ p Hamiltonian 

The Hamiltonian describing the k • p interaction 
between Fsv and Fic is well known.EJ We are only inter- 
ested in conduction-band effective masses, so we do not 
consider remote band contributions in the valence band. 
In the conduction band, both the k • p interaction with 
the topmost valence band and remote bands contribute to 
the effective mass. However, the latter terms are rather 
small, and we neglect them here. Hence reads as 
follows 




The energy reference in Eq. (|lj) is taken at the maxi- 
mum of the valence band. Note that is spherically 
symmetric, i.e., we can choose the coordinate system to 
be ex II [110], By II [112] and || [111], which is conve- 
nient for describing the L point. The only parameters we 
need to know to specify Eq. ( |l^ ) are the band gap 
and Kane's momentum matrix elementu 

= -^-{rMri^). (15) 

m 

The situation at the L point is very similar, in that 
there is only one reduced matrix element 

P"- = -*^(Lieb.|L^v). (16) 

If we neglect remote band contributions in the conduction 
band, P^ is responsible for the conduction-band trans- 
verse mass m}[. However, interactions between Lie and 
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Lsv states cannot account for the conduction-band lon- 
gitudinal mass mjj", and the longitudinal mass would be 
equal to the free-electron mass m without contributions 
from remote bands. Therefore, in the L point k • p matrix 



we have to retain the parameter G, which represents re- 
mote band contributions to mjj" . Neglecting remote band 
contributions in the valence band, the Hamiltonian ma- 
trix takes the form 



+ hlK^ + GKl iP^K, iP^Ky 
H^{K) = I -^P^K, + | . (17) 

-iP'^Ky El + 



B. Matrix elements of Vi C. Values of the matrix elements 



The potential of the ordered material can be modeled 
by dividing the lattice into two sublattices which are rich 
in Ga or In, respectively, and averaging separately over 
the two sublattices. Subtracting the averaged potential 
of the disordered material, one obtains a model for the 
ordering potential as used in Ref. |lj. In principle, there 
are two different types of matrix elements of the order- 
ing potential, those which couple F- and L-point states, 
and those which lead to interactions within T- or L-point 
states, respectively. However, the latter matrix elements 
are exactly zero, if the ordering potential is modeled with 
the above outlined separate virtual crystal approxima- 
tions over the two sublattices. Therefore, it can be ex- 
pected that these matrix elements are small, and they 
are neglected here. The nonzero matrix elements of the 
ordering potential Vi can be derived using group theory, 
and we are left with only three real reduced matrix ele- 
ments 

1^1 = (ric|Fi|Lic) (18a) 

v^35 = (r^vl^i|L3v) = (rLl^i|L|v) (18b) 

Viz = (qjFilLie). (18c) 

These equations illustrate that our generalized approach 
shares the well-known and important feature of standard 
k • p theory that by means of group theory the number 
of independent parameters can be greatly reduced. 

Combining Eqs. (|l^), ([l7|), and ( |l8| ) we end up with a 
Hamiltonian of the form 



i?^^(K) 



rL 



ffr(K) V 



(19) 



with 



/ Vii 






















^35 











(20) 



Figure ^ shows a schematic picture of the interactions 
described by the different matrix elements. 



Two limiting cases are used to determine the numer- 
ical values of the potential and momentum matrix ele- 
ments in Eq. (|l^). For Vn = V15 = V35 — the model 
describes the disordered material, and the unknown pa- 
rameters P^ , P^ and G can be fitted to the conduction- 
band effective masses at the F and L points, respectively. 
Such an analysis-using experimental data has been done 
for the F point, Ell but not for the L point. In order to 
obtain a consistent set of parameters, we deduce the ef- 
fective masses and band gaps from a band-structure cal- 
culation based on an empirical tight-binding model with 
sp^cfs* basis, nearest neighbor intecactions, and without 
spin-orbit interaction. Jancu et alE3 showed that a tight- 
binding model with such a basis is capable of accurately 
describing the valence bands and the two lowest conduc- 
tion bands in many diamond and zinc-blende-type semi- 
conductors. The tight-binding parameters we use are in- 
terpolated from the values for GaP and InP in Ref. |2|, 
with a Ga:In ratio of 51:49. In order to correctly repro- 
duce the fundamental band gap in this virtual-crystal 
approximation, we incorporate an empirical bowing fac- 
tor = l/2[(/Z'm)GaP — (^^'TO)inp] for the four (ssa)- 
type two-center integrals. Band gaps and effective masses 
from this calculation and the resulting values for P^, P^ 
and G are summarized in Table |. We use a phase con- 
vention for the wave functions, such that both P^ and 
P^ are positi\p. The values for P^ and P^ are very close 
to each other ,E3 so we use the approximation 



P^ =P^ = 8.86 eVA. 



(21) 



With nonzero potential matrix elements but K = 0, 
the Hamiltonian matrix ( p^ ) describes the zone cen- 
ter states of ordered GaInP2. These states—ha-vc, been 
studied previously, both experimentallyli3llj'oE3 and 
theoretically.EjO These studies indicate that there is 
a certain correlation between different ordering-induced 
changes of the band structure. In particular, the crystal- 
field splitting AcF, the band-gap reduction Aii^BGR, 
and the change in the transition energy Ai?r-»L for the 
ordering-induced transitiop-r3v(F5v) ric(Lic) have a 
fixed ratio for all samples :E2l 



AE'bgR „ 

C = — = 2.66, 



A 



CF 



(22a) 
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AE- 



= 0.48. 



BGR 



(22b) 



Expressing Aqf, AE'bgr, and AiJr^L as functions of 
Vii, Vi5, and V35 to second order in these matrix ele- 
ments and using ratios (p^), we obtain 



^35 



E}; 



-0.426 
-0.110 



^11 



- EF 
1 



E^-Er 



(23a) 
(23b) 



Equation (|2^) determines the relation between \Vii\ 
and IV35I, and between |Vii| and jVisj. Different degrees 
of ordering, i.e., different strengths of the ordering po- 
tential, can therefore be modeled by different values of 
Vii- The matrix element Vn itself is proportional to the 
degree of ordering 77, as defined in Ref. |lj, if the or- 
dering potential Vi{r) is described by separate virtual- 
crystal approximations over two sublattices described 
above. Note that this method does not determine the 
relative signs of the matrix elements in (pi 



D. Diagonalization 

The band-gap redjictian in highly ordered samples 
is about 150 meV.llall3'll3 This corresponds to Vn 
200 meV in our model. Thus, according to Eq. (23b), 
the potential matrix element V15, which couples Lie and 



states, is small compared to the energy difference 
between these states. We therefore use Lowdin per- 
turbation theoryna to calculate the change in energy of 
these states to second order. For the Lie state this gives 
Ec + I^i5p/-^J"; whereas for the T^^ state the energetic 
position of the level is E^^ — —\Vi5\'^/E^. Neglecting 
the mixing of wave functions, we decouple valence and 
conduction band with respect to the ordering potential 
by this procedure. 

The problem thus reduces to two two-level systems, 
which can be solved analytically, resulting in energy 
eigenvalues E^^^^^ and Ei^^^\ and expansion coefficients 
for the zone-center states in the conduction band. 



|fle(rie)) =ae|rie)+/3e|Lle), 
|fle(Lie)) =/3e|rie)-ae|Lle), 



and in the valence band. 



flv(rfv)> 



|f?v(L3v)> 



|rf,)-f/3v|L^j, 



ir 



5v/ 



a, 



ILL) 



(24a) 
(24b) 



(25a) 
(25b) 



In these equations a bar denotes states of the ordered 
material. In addition, the main contributing state of the 
zinc-blende crystal is given in parentheses. As states ( |2^ ) 
and (^H) are diagonal with respect to the ordering poten- 
tial, this removes the potential matrix elements from the 
Hamiltonian (19), but at the price of introducing new 



k • p interactions. The following four momentum matrix 
elements appear: 





h 

—i — 
m 


(fie(rie)b.|f^,(riJ) 


= -i — {T 
m 


Pl = 


n 

—i — 
m 


(fle(rie)b.|f^,(L^J) 


h 

= -i-{T 
m 




h 

~i — 
m 


(fie(rie)b.lfiv(riJ) 




p! = 


n 

~i — 
m 


(fle(Lie)b.lfiv(r^J) 





(26a) 
(26b) 
(26c) 
(26d) 



where we have already used relation (^ij). The momentum matrix elements (|2^) define a standard k ■ p problem of 
the form of Eq. (0) for ordered GaInP2, 



P / i/r(r)(K) i/r(r)r(L)(K) 



(27) 



with 



J?^(L)(K) = 



E. 



(2) 



-iP.^K, 

~^PlKy 

-~iP}k, 

-iP^K^ 
-iPtKy 



E. 



(1) 



2m 















EV + 






El' + €lK' J 



E, 



iPtKy 



(2) 



5 



and 



ffr(r)f(L)(j^) ^ 



/ 





-iP^K.. 
























V 


-iP^K, 





I 



r 



Without the approximation of Eq. (^Tj) the form of 
Hamihonian ( p7| ) would correspond to the general case 
of a crystal with C^^ symmetry. The momentum ma- 
trix elements P^ and determine the effective masses 
of ric(ric) and ric(Lic) perpendicular to the ordering 
direction. The momentum matrix elements pj' and p]' 
together with G determine the effective masses parallel 
to the ordering direction. A schematic picture for the in- 
teractions described by the momentum matrix elements 
is shown in Fig. 0. 



IV. RESULTS AND DISCUSSION 

Having set up our model, we can first calculate \Vi^\ 
and iVijsl for different values of Vn, and then derive the 
new band-edge energies and the expansion coefficients 
in Eqs. ( |2^ ) and ([25|). The expansion coefficients deter- 
mine the momentum matrix elements (p^), which, to- 
gether with the new band-edge energies, yield the effec- 
tive masses. The results for the momentum matrix ele- 
ments and effective masses are plotted for a range of | Vli | 
up to 0.35 eV. This value results in a band gap reduction 
Ai^BGR of about 430 meV, which is the theoretical value 
for the perfectly ordered CuPt structure.EJ 

Up to now we have not considered the different possi- 
bilities for the relative signs of the potential matrix ele- 
ments. The sign of Vx^ does not matter since only |Vi5|^ 
enters into a second order perturbation theory correc- 
tion. Hence only the relative sign of Vu and V35, that is 
a = sgn(Vii/V35), has to be determined. We will show 
that this can be done by appropriate comparison with 
experimental results. 

The results for the squares of the four momentum ma- 
trix elements (|6|) are shown in Fig. ||(a) for ct < and 
in Fig. ^b) for cr > 0. The intensity of the optical 
transition f3v(r5v) f ic(ric) is proportional to jPf^P, 
whereas the intensity of the ordering-induced transition 
f3v(r5v) fic(Lic) is proportional to |P2"'"P- Ex- 
perimental results indicate that the latter transition is 
much weaker than the former, even for highly ordered 
samples.E^ Therefore, we can rule out the option a > 0, 
as this would result in approximately the same intensity 
for these two transitions. 

.11 



The difference between P^ and Pj 



and pj' should influence the optical anisotropy of or- 
dered GaInP2. ftliSi effect has been neglected in previ- 
ous calculations .c3E3 The mirror symmetry in Figs. 0(a) 
and||(b) with respect to a horizontal line at {P/P^Y = 



and between P^ 



0.5 is 

states 

{pI? 

not a 
hence 



due to the 
Ml) and (E5I) 



normalization of the zone center 
. The matrix elements (Pj')^ and 
do not depend on ct, since the f iv(r5y) state is 
mixture of two different zinc-blende states, and 
there are no "interference" terms in Eqs. (26c) 
and (|26d|) . In Fig. |(b) the relation {P^f « {P^f and 
(PjJ-)2 ~ (pj)2 foj. = 0.35 eV are purely accidental. 

Figure ^(a) shows the effective masses of the low- 
est conduction-band state fic(ric) state for ct < 0. 
The effective mass parallel to the ordering direction 
m|| increases with ordering, whereas the effective mass 
perpendicular to the ordering direction m±_ decreases. 
Within our model the anisotropy of the effective masses 
is (mil - mj.)/m^ = 0.489 for \Vii\ = 0.35 eV, i.e., for 
perfect ordering. This value is in good agreement with 
the results of Ref. ^ The general trend of the increasg 
in m|| and reduction of m± agrees with both theoreticalQ 
and experimental^ results. For completeness Fig. ^(b) 
shows the effective masses for ct > 0. It illustrates how 
important it is to determine ct correctly. 

The predictions of our model for the effective masses of 
the f ic(Lic) state are shown in Figs. |^(a) and^(b), again 
for CT < 0. The most striking feature is the decrease in 
the effective mass parallel to the ordering direction from 
1.7 to less than 0.4. The effective mass perpendicular 
to the ordering direction shows an increase, comparable 
in magnitude to the changes for the ric(ric) effective 
masses. To the best of our knowledge, the present work 
is the first investigation of the effective masses of this 
second lowest conduction band in the ordered material. 
For completeness Figs, ^(c) and ^(d) show the effective 
masses for ct > 0. Note that m\\ does not depend on ct 
m Figs. I and |. This can be easily understood, as these 
masses are determined by the terms proportional to G in 



Eq. (|23) and by (Pj')^ or {P2Y , respectively, which are 
independent of ct. 



V. SUMMARY AND OUTLOOK 

In conclusion, we have presented a general formalism 
which extends standard k • p theory to periodic pertur- 
bations, which are rapidly oscillating on a length scale of 
a few lattice constants. We choose a suitable complete 
and orthonormal basis that makes it possible to consider 
explicitly the interactions due to the perturbation. Our 
ansatz can be readily combined with other extensions of 
k • p theory such as for the inclusion of strain and spin- 
orbit interaction, thereby retaining the simple analytic 
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formulas of k • p theory. When the period of the per- 
turbation increases, more K. points in the Brillouin zone 
have to be considered in our ansatz, increasing the num- 
ber of potential matrix elements . However, the 
number of independent parameters can be significantly 
reduced using symmetry arguments, as illustrated in our 
discussion of ordered GaInP2. If it is desirable to de- 
crease the number of parameters further, one can in a 
perturbative way restrict the calculation to a subset {K} 
containing only extremal points of the energy dispersion, 
which usually are most important. Alternatively, one 
can calculate the potential matrix elements according to 
their microscopic definition [Eq. (|l|)] using wave fuctions 
from, e.g., a pseudopotential calculation for the unper- 
turbed system. 

As an example, we calculate the effective masses in the 
lowest two conduction bands of spontaneously ordered 
GaInP2 as a function of the degree of ordering. For the 
lowest conduction band we find qualitatively good agree- 
ment between our results, first-principle calculations and 
experimental data. We also find the momentum matrix 
element between conduction- and valence-band states to 
be anisotropic, which influences the optical anisotropy of 
ordered GaInP2. Although we have calculated the curva- 
tures of the conduction bands only, our approach can also 
be applied to the valence band. To do this, a consistent 
set of band parameters is required and spin-orbit interac- 
tion should be taken into account. We expect, e.g., that 
the different signs of the curvature of the valence band 
parallel to the ordering direction at T and L points cause 
an increase in the heavy-hole mass parallel to the order- 
ing direction. Besides GaInP2, which we have treated 
here, natural short-period superlattices occur in many 
different semiconductor alloys (cf., e.g., Ref ^4|), and our 
method is well suited to describe these systems. 
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FIG. 1. Schematic picture of the interactions described by 
the momentum and potential matrix elements in Eq. 




FIG. 2. Schematic picture of the interactions described by 
the momentum matrix elements after diagonalization with re- 
spect to the ordering potential [Eq. (^^]. For clarity the pic- 
ture is split into (a) K perpendicular and (b) K parallel to 
the ordering direction. 
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FIG. 3. Square of momentum matrix elements [Eq. (p6[)] 
for (a) a = sgn(Vii/F35) < and (b) cr > 0. 
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FIG. 4. Effective masses of the lowest conduction-band 
state ric(ric) for (a) cr = sgn(Vii/V35) < and (b) cr > 0. 
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FIG. 5. Effective masses of the second-lowest conduc- 
tion-band state fic (Lie) for (a) and (b) cr = sgn(Vii/p36) < 0, 
and (c) and (d) cr > 0. 

TABLE L Energies, effective masses, and momentum ma- 
trix elements or remote bands contribution, respectively, for 
Tie, Lie, and Lsv The energy zero is the valence-band max- 
imum Fsv. 



State 


Energy 


Effective mass 


Matrix element 




(cV) 


(m) 




Lie 


E'^ = 2.024 


m'' = 0.0899 


P' = 8.83 eVA 


Lie 


E^ = 2.250 


mi = 0.1349 


= 8.88 eVA 






mj| = 1.699 


G = -1.57 eVA^ 


Lsv 


E^ = -0.978 
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